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Abstract — A  completely  new  type  of  solid-state  laser  device  is 
presented  that  offers  the  potential  for  achieving  significantly  in¬ 
creased  levels  of  terahertz  (THz)  frequency  output  power  at  rela¬ 
tively  high  operating  temperatures.  Specifically,  a  double-barrier 
GaSb/InAs/GaSb  heterostructure  device  concept  is  introduced  that 
simultaneously  leverages  resonant  electron  injection  and  interband 
tunneling  electron  depletion  to  realize  electron-population  inver¬ 
sion,  while  at  the  same  time  mitigating  the  scattering  effects  that 
degrade  the  lasing  process.  Here,  the  main  innovations  are  the  abil¬ 
ity  to  spatially  separate  the  upper  and  lower  electron  populations 
using  the  quantum  confinement  of  the  double-barrier  conduction 
band  well  and  the  valence-band  (VB)  well  (i.e.,  of  the  second  bar¬ 
rier),  respectively,  and  the  depopulation  of  the  VB  well  by  heavy 
hole  interband  tunneling.  A  theoretical  analysis  of  the  radiative 
and  nonradiative  transition  rates  based  upon  a  multiband  Kane 
model  formalism  is  used  to  confirm  the  large  available  optical  gain 
and  to  estimate  the  lasing  output  power  at  very  long  wavelengths. 
Therefore,  this  study  establishes  the  initial  foundation  for  a  solid- 
state  THz  laser  that  can  provide  significant  levels  of  output  power 
below  1  THz.  Furthermore,  the  inherently  high  spectral  purity  and 
natural  tunability  offered  by  this  novel  laser  technology  will  be  in¬ 
strumental  in  future  spectroscopic  analysis  of  nanoscale  biological 
and/or  organic  systems  that  are  well  known  to  possess  unique  spec¬ 
tral  signatures  at  very  long  wavelengths. 

Index  Terms — Conduction  band,  heavy  hole,  heterostructure, 
InAs/GaSb,  interband,  lasing,  nonradiative  transitions. 


I.  Introduction 

THE  key  technological  property  of  broken-gap  InAs/GaSb 
heterostructure  from  the  6.1  A  family  is  that  the  valence- 
band  (VB)  bottom  of  GaSb  is  0.15  eV  above  the  conduction- 
band  (CB)  bottom  of  InAs  [1],  Hence,  within  the  adjacent  InAs 
conduction-band  (CB)  well  and  the  GaSb  heavy-hole  (HH)  VB 
well  systems,  the  individual  electron  and  hole  energy-levels  may 
be  positioned  close  enough  to  allow  for  very  long  wavelength 
recombination  while  the  same  time  the  overlaps  of  individu¬ 
ally  quantum  confined  wave  functions  are  small.  It  is  useful  to 
note  that  this  set  of  quasi-bound  systems  allows  for  a  finite  in- 
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terband  (i.e.,  conduction  to  valence)  recombination  process  to 
exist  and  that  the  small  overlap  will  significant  reduce  important 
scattering  effects  (e.g.,  phonon  and  Auger)  that  typical  degrade 
lasing.  Furthermore,  as  will  be  shown,  it  is  possible  to  utilize 
interband  tunneling  to  rapidly  drain  the  electrons  from  the  VB 
well  and  to  achieve  significant  population  inversions  and  optical 
gain.  The  fact  that  alloying  with  aluminum  or  indium  may  mod¬ 
ify  band-edge  alignments  adds  additional  flexibility  for  design 
of  the  device’s  operational  frequency  and  lasing  performance. 
For  these  general  reasons,  the  InAs/GaSb  heterostructure  is  of 
great  interest  for  applications  in  long  wavelength  light  emit¬ 
ting  devices  and  detectors.  To  date,  these  investigations  have 
been  limited  to  the  study  of  radiation  sources  in  the  infrared 
regime,  with  some  examples  being:  Type-II  intersubband  lasers 
as  described  in  [2]  and  [3];  a  Type-II  cascade  lasing  device 
as  given  in  [4];  and,  a  W-shape  Type-II  midinfrared  laser  dis¬ 
cussed  in  [5].  This  paper  extends  a  more  recent  development  of 
a  unipolar  double-barrier  staggered  AlGaSb/InAs/AlGaSb  in¬ 
terband  tunneling  diode  (I-RTD)  for  midinfrared  lasing  [6],  [7]. 
The  basic  concept  is  to  utilize  the  electron  injection  produced 
by  the  resonant  tunneling  through  the  CB-well  structure  (i.e., 
formed  by  the  double  AlGaSb  barriers  and  InAs  well)  as  the 
source  of  the  lasing,  and  the  HH  population  generated  within 
the  VB  well  (i.e.,  formed  by  the  right  AlAsSb  barrier)  by  in¬ 
terband  Zener  tunneling  as  the  recombination  sink.  The  basic 
utility  of  the  broken-gap  I-RTD  lasing  process  for  realizing  far- 
infrared  to  THz  radiation  was  recently  summarized  by  our  group 
in  [8],  and  the  goal  of  this  study  is  to  assess  a  specific  InAs/GaSb 
broken-gap  heterostructure  laser  design  that  has  been  tailored 
for  THz-frequency  operation. 

The  basic  energy-band  structure  of  the  broken-gap  I-RTD 
laser  device  considered  in  this  study  is  illustrated  in  Fig.  1(a). 
As  illustrated  in  Fig.  1  (b),  the  left-hand  side  of  the  device  will 
be  implemented  with  a  heavily  doped  n+  InAs  layer  on  the  left 
which  forms  the  diode  emitter  and  a  combined  undoped  (spacer) 
and  n+  InAs  layer  on  the  right-hand  side,  which  forms  the  diode 
collector.  Electrons  will  then  be  injected  into  the  higher  energy 
level  E2  that  is  quasi-confined  in  the  InAs  CB-well  region  when 
the  structure  is  subjected  to  a  forward-biasing  external  voltage, 
and  this  will  be  accompanied  by  a  depopulation  of  electrons 
from  the  lower  energy  level  E\  that  is  quasi-confined  in  the  right 
GaSb  VB  well  region  by  HH  interband  tunneling  of  electrons 
into  the  undoped  spacer  region. 

The  carrier  inversion  for  the  two-energy  system  defined 
in  Fig.  2  may  be  estimated  using  the  equation  [9],  AN  = 
ViJ(T2i  ~  ri)/e>  where  rji  is  the  quantum  efficiency,  J  is  the 
CB  resonant  current,  T21  is  the  relaxation  time  constant  of 
energy  level  2,  and  iq  is  the  time  constant  for  removal  of 
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Fig.  1 .  Broken-gap  device  structure  for  THz  lasing,  (a)  Without  bias,  (b)  Under 
bias. 


carrier  from  energy  level  1.  The  relaxation  time  is  written 
as  1/ T21  =  1/tsp  +  1/Tap  +  1/Top  +  l/TAug,  where  rsp  is  the 
spontaneous  recombination  time  constant,  rap  is  the  acoustic 
phonon  scattering  time  constant,  rop  is  the  optical  phonon  scat¬ 
tering  time  constant,  and  TAUg  is  the  Auger  recombination  time 
constant.  The  key  to  establishing  a  useful  population  inversion 
is  that  charging  rate  (1/ri)  of  the  E\  by  HH  tunneling  can  be 
made  much  faster  than  the  net  relaxation  rate  (I/T21).  This  fact 
will  be  substantiated  by  the  modeling  and  simulation  work  that 
is  presented  shortly. 

Section  II-A,  II-B,  and  II-C  will  present  physical  models  that 
are  accurate  for  addressing  the  most  important  issues  related 
to  the  broken-gap  I-RTD  laser  device;  namely,  the  multiband 
transport  models,  transition  probabilities,  and  the  self-consistent 
charging  effects  which  can  be  combined  together  to  estimate 
all  the  flux  processes  and  the  overall  depopulation  dynamics. 
Section  II-D  provides  an  assessment  of  a  specific  broken-gap 
I-RTD  laser  design  that  targets  THz  frequency  operation.  The 
main  conclusions  of  this  study  are:  1 )  CB  energy  level  Ei  and 
heavy  hole  level  E\  are  amenable  to  the  production  of  long 
wave  radiation;  2)  the  HH  interband  tunneling  is  fast  enough  for 
electron  depopulation;  3)  nonradiative  process  rates  that  degrade 
performance  can  be  suppressed;  and  4)  substantial  optical  gain 
can  be  obtained. 

II.  Modeling  and  Simulation  Results 

The  study  of  transport  dynamics  in  this  broken-gap  and 
narrow-band  structure  must  be  based  upon  a  six-band  Kane 


Injection  by  CB  resonant  tunneling 
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simulated  emission 


■ 


AA. 


N on-radiative  process 


E, 


Depopulation  by  HH  interband  tunneling 


Fig.  2.  Two-energy  level  model  for  electron’s  injection  and  removal.  The 
energy  levels  in  emitter’s  and  collector’s  reservoir  are  not  sketched. 


model  formulism  of  the  type  previously  developed  in  [10] 
and  [11].  This  model  allows  for  the  treatment  of  the  strong 
multiband  coupling  effects  that  originate  from  the  narrow 
bandgaps  and  band  edge  misalignments.  A  complete  descrip¬ 
tion  of  the  CB  electron  motion  requires  a  set  of  envelope  func¬ 
tions  t/08^  =  |cicaB* ,  j,  where  ai88'1  is  the  conduc¬ 

tion  band  (CB),  is  the  light-hole  (LH)  band,  and  a'^  is 
the  heavy-hole  (HH)  band  components.  Here,  the  subscript  a 
denotes  the  two  possible  spin  state,  i.e.,  ±.  Similarly,  the  HH 
electron  dynamics  in  the  VB  requires  another  set  of  envelope 
functions  denoted  as  t/™-*  =  |ai8H') ,  ,  a j. 

The  net  CB  transmission  of  electrons  from  emitter  to  collec¬ 
tor  arises  through  contributions  from  two  distinct  channels.  The 
first  channel  is  from  resonant  tunneling  through  the  CB-well 
that  is  defined  by  the  GaSb/InAs/GaSb  double-barrier  struc¬ 
ture.  This  electron  flux  component  is  illustrated  by  the  upper 
horizontal  arrow  in  Fig.  1(b).  The  second  channel  is  from  in¬ 
terband  tunneling  of  electrons  that  are  allowed  to  leave  the  HH 
level  within  the  right  VB  well  and  that  enter  into  the  spacer- 
collector  region.  This  electron  flux  component  is  illustrated  by 
the  lower  horizontal  arrow  on  the  right  side  of  HH-state  E1  as 
also  shown  in  Fig.  1(b).  Here,  it  is  valuable  to  note  that  these  VB 
well  electrons  can  be  sourced  from  interband  tunneling  from  the 
emitter  [i.e.,  see  lower  horizontal  arrow  on  the  left-hand  side  of 
HH-state  E\  in  Fig.  1(b)]  and  from  electron  recombination  that 
occurs  from  f?2  to  E\  transitions.  However,  as  will  be  demon¬ 
strated  later,  the  contribution  from  interband  tunneling  from  the 
emitter  is  insignificant  when  forward  biasing  is  applied.  This 
means  that  HHs  will  accumulate  in  the  right  VB  well,  and  then 
collectively  serve  as  a  localized  recombination  center  for  CB 
electrons  that  are  injected  into  the  CB-well.  This  recombination 
will  arise  from  both  radiative  or  nonradiative  processes. 


A.  Transport  Equations 

The  Kane  model  from  [10]  and  [11]  is  useful  for  describing 
the  wavefunctions  through  the  various  regions  of  the  broken- 
gap  structure  illustrated  in  Fig.  1.  The  CB  electrons  in  the  InAs 
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regions  are  governed  by  the  coupled  differential  equations 


1. 
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where  Egw  is  the  band  gap  of  InAs,  V„  is  the  bias  voltage,  kt  = 
P  ky  is  the  in-plane  momentum,  e  is  the  charge  of  electron, 
and  P  is  Kane  interband  parameters.  Here,  E  is  referenced  to 
the  CB  bottom  in  the  emitter.  The  coupled  differential  equation 
system  for  CB  electrons  in  the  GaSb  regions  is 
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Here,  77  is  the  reflection  coefficient  and  N)  is  the  normalization 
constant  [10].  In  both  (8)  and  (9),  the  first  term  represents  the 
incident  wave  beginning  from  the  emitter  contact  while  the  sec¬ 
ond  term  accounts  for  the  reflection  by  potential  barriers,  also, 
kRz  is  the  wave  vector  of  propagation.  Therefore,  the  system 
defined  by  (1),  (2),  (8),  and  (9)  leads  to  a  dispersion  law  of 

E(E  P  Egw)  =  (*L  +  k2t)  (10) 

where  mw  =  0.023mo  is  the  electron’s  effective  mass  in  InAs 
layers.  Furthermore,  if  the  resonant  energy  level  in  the  CB-well 
is  referenced  to  the  VB  edge  at  the  collector  contact  (this  will 
be  needed  later)  then  it  can  be  expressed  as 

Ei  =  E  P  Egw  P  eVa.  (11) 

The  interband  tunneling  can  only  through  CB-LH-HH  multi¬ 
band  coupling  effects,  and  it  is  possible  to  deduce  that  it  will  be 
dominated  by  the  HH  confined  state  within  the  right  VB  well 
(i.e.,  HH  state  will  be  closet  to  the  VB  edge).  Hence,  it  is  accept¬ 
able  to  assume  that  a{HH*  and  ai™*  are  an  order  of  magnitude  less 
than  a\  ,  which  leads  to  the  following  systems  of  differential 
equations  describing  the  wavefunctions  in  the  right  VB  well 
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where  Egb  is  the  band  gap  of  GaSb;  and  Ac  is  the  conduction 
band  offset.  Due  to  the  fact  that  the  LH  and  HH  contributions 
to  the  CB  transport  should  be  small,  the  Kane  model  directly 
reveals  that  the  HH  component  may  be  estimated  from  the  linear 
equation 
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For  resonant  transmission  in  the  CB,  the  wavefunction  in  the 
collector  on  the  right  ( R )  is  only  an  outgoing  wave,  which  leads 
to  the  relations 


(cb)  _ +  exp  {ikRzz) 
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where  tRc  is  the  transmission  coefficient,  Nr  is  the  normaliza¬ 
tion  constant,  and  the  outgoing  wave  vector  kR  can  be  obtained 
from  the  dispersion  relation  of  energy  [10].  The  wavefunction 
in  the  emitter  on  the  left  ( L )  must  be  the  sum  of  incoming  and 
reflected  waves,  which  leads  to 


a 


(CB) 

Leo 


1 


[exp  ( ikLzz )  +  rL  exp  {-ikLzz)\ 


(8) 


=  -(Egb  PUP  Av  -  Eh)a\ W  ±  ^Pkta{™]  (12) 

\/IP  [j~z  ±  H  =  {Egb  +  U  +  A"-  <13) 

&h^-lEh-Uiz)-A^=0  (14) 

where  A„  is  the  VB  offset,  rrih  =  0.40mo  is  the  effective  mass  of 
heavy  hole,  and  Eb  is  the  quasi-bound  heavy-hole  level.  All  the 
energies  are  referenced  to  the  CB  bottom  at  the  collector  contact. 
In  the  InAs  regions,  the  differential  equations  are  obtained  by 
replacing  the  band  gap  Egb  with  Egw ,  and  letting  Av  =  0.  The 
localized  heavy-hole  wavefunctions  at  the  boundaries  of  the 
device  contacts  can  be  written  as 

aiH£TH)  =  bh.L  exp  ( KLhz )  (15) 

=  bhR  exp  (- KRhz )  (16) 

where  nLh  =  yj2 mh{Eh  -  eVa)/Ti  and  nRh  =  y/2mhEh/h  at 

the  emitter  and  collector  contacts,  respectively.  The  quasi-bound 
Eb  can  be  estimated  from  ( 14).  The  solutions  of  and  a'™1 
in  the  emitter  region  must  be  written  as 

affi  =  BcLa  exp  (ikLz)  +  bcLa  exp  (klhz)  (17) 


(CB) 

’Ll 


1 

( E  P  Egw)^/NR 


exp  ( ikLzz ) 


+ 


-kLz 


rL  exp  (~ikLzz) 


(9) 


a\a  =  BiLa  exp  (ikLz)  P  hLa  exp  (nLhz)  (18) 

where  the  propagating  terms  allows  for  the  possibility  of  elec¬ 
trons  from  the  emitter  region  to  successfully  enter  into  the  HH 
level  of  the  right  VB  well  via  interband  tunneling.  Of  course, 
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these  electrons  contributes  no  radiation.  The  wave  vector  ki  is 

kL  =  [Eh  ~ {Egw + eK)]  {Eh  ~ eVa)  ~  \kl  (19) 
Similarly,  the  solutions  in  the  collector  region  are 

aJ.”H)  =  AcRa  exp  ( ikRz )  +  bcRa  exp  (• -KRhz )  (20) 

<4™  =  AlRa  exP  {ikRz)  +  bma  exp  (-nRhz)  (21) 

where  the  propagating  terms  represent  the  possibility  of  elec¬ 
trons  enter  into  collector  via  interband  tunneling.  The  wave 
vector  kR  is 


kR  = 


Eh{Eh  -  Egw )  - 


2  P2 


1 

n 


4  P2 


^r\ACR±(kt)\2  - 


kL 


- \BcL±(kt)\ 2 


E\  —  Eh  — 


2m/j 


state  is  1 1 CB) )  =  ai™  |  Sa }  +  a 


,(cb) 

'la 


Via 


,(CB) 


h±>  =  ^  S  exp  (—i<j>) 


± 


-lxjl\Z])  +  i^\{X  +  iY)  [) 


Z  I) 


76 


|(x-*r)  T) 


±exp(-2^)-^=|(X  +  iY)  t) 

+  exp(i4>)-^=i\(X  -  iY)  |) 


(28) 


(22) 


where  the  phase  angle  is  cf>  =  arctan  ( ky  / kx ) .  kx  and  ky  are  the 
projections  of  momentum  in  x  and  y  directions.  The  key  point 
to  make  here  is  that  these  multiband  wavefunctions  can  be  used 
to  show  that  the  emitted  light  will  be  circularly  polarized  (see 
examples  in  [12])  and  the  strength  of  the  optical  emission  is 
strongly  associated  with  the  interband  interaction  parameter  P. 
By  using  these  same  wavefunctions  in  the  matrix  element  for 
optical  transition  leads  to  a  time  constant  for  the  spontaneous 
transition  rate  of 


The  rate  of  HH  electron  tunneling  out  of  (and  into)  the  VB  well 
can  be  evaluated  from 


1 


4m,  „  e~  P2nr 


n4c3 


(23) 


Eh  y  Eh  —  eVa 

where  kR  (fc/ )  is  outgoing  (incoming)  wave  vector  for  the  HHs 
and  \  AcR±{kt)\2  (\BCL±(kt)\2)  is  the  interband  tunneling  prob¬ 
ability  of  the  outgoing  (incoming)  HHs.  Hence,  (23)  is  the  cen¬ 
tral  result  for  determining  the  depopulation  dynamics.  The  net 
variation  of  charge  in  the  right  VB  well  due  to  outgoing  and 
incoming  electron  tunneling  is  now  definable  from 

dQ  I"  ktdkt  1 

dt  6  J  2n  n 

where  the  last  equation  represents  a  consistent  detailed  balance- 
type  summation  over  all  the  available  HH  states.  HH  energy  is 
defined  from 


dui(w/kLo)\{  ajfH)|acCB)  )|2/c(l  -  fv) 


1  + 


27 d-(kl0+kl) 


(29) 


mw  Ea 


where  nr  is  diffraction  index,  Till  is  photon  energy,  c  is  the 
speed  of  light,  and  kRo  is  root  of  equation  Ei  —  E\  =  Till.  fc 
is  Fermi  functions  for  the  initial  electron  state  (i.e.,  equal  to 
that  in  the  emitter);  fv  is  Fermi  function  for  final  hole  state; 


(24) 


Efc  =  Er 


■E, 


gw 


(25) 


One  should  note  that  the  energy  levels  given  in  both  (11)  and 
(25)  include  the  energy  contributions  from  the  quantum  confine¬ 
ment  (i.e.,  the  --direction)  and  from  the  in-plane  direction.  This 
is  actually  very  important  as  this  interband  laser  will  derive  ra¬ 
diation  from  many  equivalent  energy  differences  that  arise  from 
these  nonlinear  dispersion  laws  and  this  lead  to  increased  optical 
gain  and  bandwidth  as  compared  to  an  intraband  laser. 

B.  Transition  Processes 

1 )  Optical  Recombination:  The  matrix  element  for  optical 
transition  is  (f\~t  ■  p|i)  between  initial  state  |i)  and  final  state 
|/),  where  T?  is  the  polarization  direction  of  the  emitted  light 
and  p  is  the  momentum  operator  for  the  electron.  The  initial 


i 'ha),  the  final  T 


eVa,  the  Fermi  energy  at  the  emitter,  kR 
is  Boltzmann  constant,  and  T  is  temperature.  The  major  nonra- 
diative  processes  include  acoustic/optical  scattering  and  Auger 
recombination  and  they  will  be  considered  in  the  following  two 
sections. 

2)  Phonon  Scattering:  The  initial  states  consist  of  CB  elec¬ 
trons  are  undergoing  resonant  tunneling  and  therefore  possess 
a  momentum  {kx,  ky ,  fc^},  where  kx  is  the  x-direction  wave 
vector  and  ky  is  the  y-direction  wave  vector.  Conversely,  the 
final  states  are  confined  HHs  and  only  possess  in-plane  mo¬ 
mentum  {kx',kyi}  where  kx'  is  the  x-direction  wave  vector, 
and  ky'  is  the  y-direction  wave  vector,  q  is  the  wave  vec¬ 
tor  of  phonon  with  projected  spatial  components  {qXl  qy:  qz}. 
q  =  jq|  is  the  magnitude  of  phonon  wave  vector.  For  phonon 
emission,  the  in-plane  momentum  satisfies  qx  =  kx  —  kx'  and 
qy  =  ky  —  ky'.  Hence,  the  magnitude  is  evaluated  from  q  = 
\J (kx  —  kx')2  +  (ky  —  ky')2  +  <?2.  The  acoustic  and  optical 
phonon  scattering  rates  for  a  specified  final  state  are  written 
as  [13] 

16t P’Tips2  J  d3fcrf9z|%)|2(/c  -  fv) 


1 


state  is  |T(HH))  =  a{™’\vha),  where  Sa),  \vi„),  and  \vha)  are 

the  Bloch  base  functions 

\S±)  =  -j=[-iexp(-i(/>)\S  T)  ±  \S  |)]  (26) 


exp 


/  TlUJq  \ 

[k^rj  ~ 


1  -1 


+  1  }  5{E2  -E1-  hujq)  (30) 


and 
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(27) 


47T2 


exp 


£o 


TlU/q 

ErT 


- 1 


+  l}6(E2-E1-huq)  (31) 
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where  D  is  deformation  potential,  Tiu>q  is  acoustic  phonon  en¬ 
ergy,  s  is  velocity  of  sound,  p  is  the  mass  density  in  the  semicon¬ 
ductor,  Tiluq  is  the  optical  phonon  energy,  s^.  is  high-frequency 
dielectric  constant,  eo  is  static  dielectric  constant.  For  the  acous¬ 
tic  phonon,  there  is  dispersion  relationship  uq  =  sq.  The  matrix 
element  is  written  as 

I(Qz)  =  J  aT *  exp  (iqz z)a{™] dz  (32) 


dynamics  to  estimate  the  total  optical  gain  as  [7] 


where  kto  need  be  evaluated  from  E2  E\  =  huj,  where  E2  is 
defined  in  (1 1)  and  £j  is  defined  in  (25). 


where  the  wavefunctions  are  calculated  from  multiband  Kane 
models  defined  in  Section  11-A. 

3)  CHHH  Auger  Recombination:  A  CB  electron  may  re¬ 
combine  without  emitting  radiation  by  releasing  kinetic  energy 
to  a  second  HH  in  the  proper  amount,  i.e.,  E2  —  E\.  In  this 
so-called  “CHHH”  process,  the  CB  electron  (i.e.,  the  “C”  in 
“CHHH”)  initially  has  the  energy  Ee y  =  eVa  +  Egm  +  E(qy) 
and  momentum  q^,  in  which  E(qi>)  is  defined  the  same  as  in 
(10).  The  recombining  HH  (i.e.,  the  first  “H”  in  “CHHH”)  ini¬ 
tially  has  the  energy  Eh i  =  Eh  —  (h2ql /2mh)  and  momentum 
q( .  During  the  Auger  process,  the  recombining  HH  will  give 
up  energy  to  a  second  HH  (i.e.,  the  second  “H”  in  “CHHH”) 
which  initially  has  the  energy  Ehy  =  Eh  —  {Ti2  q\, /2rrih)  and 
momentum  q2/.  After  the  Auger  process,  the  second  hole 
(i.e.,  the  third  “H”  in  “CHHH”)  will  have  the  energy  Eh2  = 
Eh  —  (Ti2q2/2mh)  and  momentum  q2.  The  projected  compo¬ 
nents  of  q  |  -  are  labeled  as  qyx,  qyy,  and  qyz.  Similar  notations 
are  used  for  qi,  q2,  and  q2-.  The  resulting  transitions  rate  for 
CHHH  Auger  recombination  can  be  estimated  from  [14] 


2  mw  e4 


dqvxdqi  lydq2xd'q2y 


TAug  Tl  S^TT^1 

<4? 1 1  exP  (%  Z1 ) \ahA] )  (<4 I  exp  (~kz  Zz)  |<4™} ) 


dqz 


9?  +  ~  qi’x)2  +  {qiy  -  qvy )2 


C.  Self-Consistent  Calculation 


The  resulting  potential  profile  for  the  active  region  of  the  de¬ 
vice  estimated  from  Poisson’s  equation  where  the  major  charg¬ 
ing  will  arise  from  the  accumulated  sheet  HH-charge  density 
and  leads  to  the  relation 


Ph  = 


efc?e 

2n 


(HH) 


(35) 


where  kte  is  the  maximum  wave  vector  for  HH  states  and  can 
be  derived  from 


Eh^^pfPkte^  (y  a|,HH)*ac±H)tfe^ 

n2kUy1+ y2) 

2mo 


Efv  (36) 


where  Epv  is  the  Fermi  energy  level  in  collector,  -j  and  72  are 
modified  Luttinger  parameters.  While  electron  accumulation  in 
the  CB  well  will  have  some  influence  on  the  operation  of  the 
device  (e.g.,  adjust  the  required  terminal  biasing),  the  primary 
driver  in  determining  the  overall  carrier  dynamics  and  associ¬ 
ated  transitions  is  the  VB-well  charging.  Hence,  to  reduce  the 
difficulty  of  the  numerical  simulations  for  achieving  the  first 
performance  estimates  of  the  laser  device,  the  inclusion  of  the 
CB-well  charging  has  been  reserved  for  later  investigations. 


1 

<7i'z0 


z  0 


2mw  Egm 


+ 


n  (qlx +<£,„) 

2ttiw  EgW 


-,  1/2 


f{EeV)[l- f{Ehl)} 


[1  -  f(Eh2)\  [1  -  f(Eh2,)]  (33) 


where  er  is  dielectric  constant.  qyZQ  is  found  from  the  energy 
conservation  relationship  Eey  —  Ehi  =  Epp  —  Ey 2 .  Here,  the 
screening  effect  is  not  taken  into  account.  Also,  CCCH  Auger 
recombination  (i.e.,  where  two  electrons  are  involved)  is  ignored 
because  the  electron  filling  over  the  available  CB  states  due  to 
resonant  emitter  injection  is  much  more  prolific  than  the  creation 
of  the  allowable  HH’s  in  the  VB  well  due  to  interband  tunneling. 
This  means  that  the  extent  of  electron  filling  over  CB  energies 
will  be  much  greater  than  the  extent  HH  creation  over  the  VB 
well  energies,  and  this  will  significantly  reduce  the  opportunity 
for  two-electron  scattering  (i.e.,  due  to  the  lack  of  unfilled  energy 
states  within  the  small  E2  —  E\  window)  as  compared  to  two- 
HH  scattering. 

4)  Optical  Gain:  Similar  to  (29),  Fermi’s  Golden  rule  may 
be  combined  with  the  multiband  models  for  electron  and  hole 


D.  Simulation  Results 

Numerical  simulation  studies  were  performed  to  derive  all 
the  wavefunctions  for  a  specific  design  of  the  structure  given  in 
Fig.  1.  Specifically,  a  device  with  the  following  semi-optimized 
structural  parameters  and  operational  conditions:  left  barrier 
width  of  2  nm,  well  width  of  12  nm,  right  barrier  width  of 
5  nm,  Fermi  level  in  emitter  of  0.21  eV,  Fermi  level  in  collector 
of  0.099  eV,  bias  voltage  of  0.093  V,  and  operating  temperature 
of  280  K.  The  resulting  CB  resonant  level  E2  and  HH  level 
E\  are  plotted  together  in  Fig.  3.  The  figure  clearly  shows  the 
transmission  coefficient  is  a  very  good  fit  to  the  Breit-Wigner 
formula  |f^c|2  =  TrT2 /[(E  —  E2f)2  +  F2] ,  where  E2q  is  the 
position  of  peak,  Tp  is  the  height  of  peak,  and  F  is  the  half 
width  of  half  maximum.  In  the  example  shown  in  Fig.  3,  Tj> 
is  0.61,  E2q  is  0.0386  eV,  and  T  is  0.0014  eV.  The  resonant 
tunneling  time  is  0.47  ps  from  rtr  ~  h/T.  The  distance  between 
the  HH  level  Eh  and  E2q  at  kt  =  0.05  is  5.2  meV  (238  pm, 
and  1.26  THz).  The  amplitude  Tp  and  the  position  of  E2q  are 
found  strongly  dependent  on  kt.  As  kt  increases,  Tp  decreases 
and  E2 0  increases. 
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Fig.  3.  Conduction  transmission  coefficient  at  kt  =  0.05/nm.  (a)  Calculated 
results,  (b)  Fitting  by  Breit-Wigner’s  function,  and  (c)  Position  of  heavy-hole 
level  Eh  relative  the  CB  edge  in  the  emitter. 


Fig.  6.  Transition  rates. 


Fig.  7.  Total  gain  at  photon  energy.  The  vertical  axis  is  in  log  scale. 


Fig.  4.  Heavy-hole  interband  tunneling  coefficients. 


Fig.  5.  Heavy-hole  interband  tunneling  time. 

The  tunneling  probabilities  \ACR±(kt)\2  and  |-Bci±(A)l“ 
given  in  Fig.  4  demonstrate  that  HH  interband  tunneling  is 
dependent  on  in-plane  momentum.  An  important  observation 
from  Figs.  3  and  4  is  that  the  interband  tunneling  probability 
is  significantly  less  than  the  CB  resonant-tunneling  probabil¬ 
ity.  The  huge  differences  suggest  the  existence  of  two  distinct 
electron  transport  pathes:  CB  resonant  transport  and  HH  inter¬ 
band  tunneling.  The  small  interband  tunneling  probability  also 
justifies  that  the  assumptions  made  in  (12)— (14).  Furthermore, 
the  plot  in  Fig.  5  indicates  right-hand  side  interband  tunneling 
dominates  over  left-side  interband  tunneling  (i.e.,  verifies  that 


(23)  is  positive)  when  the  active  device  region  is  tilted  under 
bias  voltage,  so,  HH  hole  accumulation  occurs.  In  addition,  the 
interband  tunneling  rate  is  always  significant  and  only  varying 
from  on  the  order  of  30  ps  at  kt  =  0.01/ nm  to  as  short  as  0.4  ps 
at  kt  =  0.1/nm. 

All  the  relevant  scattering  processes  and  the  optical  recombi¬ 
nation  rate  are  plotted  in  Fig.  6  as  a  function  of  the  particular 
final  state  (with  momentum  kt).  Here,  one  can  observe  that  the 
optical  phonon  scattering  lifetime  is  on  the  order  of  nanoseconds 
and  the  acoustic  phonon  scattering  lifetime  is  on  the  order  of 
hundreds  of  nanoseconds.  Also,  the  Auger  recombination  life¬ 
time  is  on  the  order  of  nanoseconds.  These  results  sufficiently 
demonstrate  that  HH-interband  tunneling  is  much  faster  than 
nonradiative  processes,  which  is  the  key  for  the  setup  of  stimu¬ 
lated  emission.  Furthermore,  since  the  optical  recombination  is 
on  the  order  of  20-100  ns,  the  interband  tunneling  is  sufficient 
to  maintain  a  population  inversion. 

The  total  gain  G  is  plotted  in  Fig.  7.  G  is  the  integral  of  optical 
gain  g(z)  over  the  active  region  which  has  the  dimension  of 
~20  nm.  Therefore,  the  optical  gain  g(z)  of  this  laser  concept 
is  significant  (~  103  cm  /)  in  THz  frequency  range  and  shows 
promise  for  providing  significant  levels  of  output  radiation  of 
high  spectral  purity  below  1  THz,  where  other  proven  laser 
concepts  (e.g.,  quantum  cascade)  encounter  difficulties.  Hence, 
these  preliminary  research  investigations  offer  a  new  avenue 
for  the  potential  realization  of  a  completely  new  type  of  THz 
solid-state  laser.  This  envisioned  THz  laser  has  relevance  to 
nanoscale  science  and  technology  because  tools  of  this  type 
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will  be  needed  for  the  comprehensive  analysis  of  nanoscale 
biological  and/or  organic  systems  that  are  well  known  to  possess 
unique  spectral  signatures  at  very  long  wavelengths  [15]. 

For  practical  laser  operation,  the  CB-resonant  tunneling  in¬ 
jection  should  be  minimized  regarding  issue  of  thermal  dissipa¬ 
tion.  The  operation  current  can  be  controlled  in  a  tolerable  range 
without  significantly  decreasing  optical  gain.  For  example,  one 
can  reduce  the  current  by  reducing  the  doping  density  (or  the 
Fermi  level)  in  emitter.  Another  possible  approach  to  reduce  in¬ 
jection  current  is  to  properly  increase  the  thickness  of  left  GaSb 
barrier. 

III.  Conclusion 

This  paper  has  provided  physical  models  and  numerical  sim¬ 
ulation  results  which  support  the  premise  that  population  in¬ 
versions  can  be  established  between  a  CB-resonant  tunneling- 
injection  channel  and  a  HH-interband  depopulation  channel  that 
occur  within  InAs/GaSb  broken-gap  heterostructures  and  allow 
for  radiative  emissions  at  terahertz  (THz)  frequencies.  Further¬ 
more,  a  physical  analysis  based  upon  the  application  of  the 
multiband  Kane  model  for  a  particular  InAs/GaSb  broken-gap 
heterostructure  laser  design  with  semi-optimized  structural  and 
operational  parameters  suggests  that  significant  optical  gain  can 
be  obtained  below  1  THz  at  high  operating  temperatures.  There¬ 
fore,  this  novel  interband  resonant-tunneling  device  (I-RTD)- 
based  laser  concept  has  important  relevance  to  the  very  long 
wavelength  spectroscopic  characterization  of  nanoscale  biolog¬ 
ical  and/or  organic  systems. 
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